class particula
{
   private:
	double xxx;
      
   public:
	   double pos[3];
	   double posOld[3];
	   double vel[3];
	   double force[3];
	   double forceR[3];
	   double forceK[3];
	   double drg[3];
	   double drv[3];
	   double rstart[3];

	   int nvecinos;
	   int lstVecinos[1000];

	   particula()
	   {

	   }

	   particula(double rr[3])
	   {
		   for (int i=0;i<3;i++)
		   {
			   pos[i]=rr[i];
			   posOld[i]=rr[i];
			   rstart[i]=rr[i];
		   }

	   }

//	   void addVecino(int ind)
//	   {
//		   lstVecinos[nvecinos]=ind;
//		   nvecinos++;
//	   }
	   void setPosition(double rr[3])
	   {
		   for (int i=0;i<3;i++)
		   {
			   pos[i]=rr[i];
			   posOld[i]=rr[i];
			   rstart[i]=rr[i];
		   }
		   xxx=rr[0];
	   }

	   void setVelocities(double vv[3])
	   {
		   for (int i=0;i<3;i++)
			   vel[i]=vv[i];
	   }

	   void setForces(double ff[3])
	   {
		   for(int i=0;i<3;i++)
		   {
			   force[i]=ff[i];
		   }
	   }

	  void boundary(double lbox)
	  {
		  for (int i=0;i<3;i++)
		  {
			  pos[i]-=floor(pos[i]/lbox)*lbox;
		  }
	  }

};

